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ABSTRACT 

We present a study of the long-term evolution of the cloud of aerosols produced in the 
atmosphere of Jupiter by the impact of an object on 19 July 2009 [Sanchez-Lavega et 
al.. 2010. The Impact of a Large Object on Jupiter in July 2009, ApJL, 715 , L155-1591- 
The work is based on images obtained during 5 months from the impact to 31 December 
2009 taken in visible continuum wavelengths and from 20 July 2009 to 2$ May 2010 
taken in near-infrared deep hydrogen-methane absorption bands at 2. 1-2.3 microns. The 
impact cloud expanded zonally from ~ 5000 km (July 19|fo 225,000 km (29 October, 
about 180° in longitude), remaining meridionally localized within a latitude band from 
53.5°S to 61.5°S planetographic latitude. During the first two months after its formation 
the site showed heterogeneous structure with 500-1 000 km sized embedded spots. 
Later the reflectivity of the debris field? became more homogeneous due to clump 
mergers. The cloud was mainlfliispersed in longitude by the dominant zonal winds and 
their meridional shear, during the initial stages, localized motions may have been 
induced by thermal per®bation caused by the impact’s energy deposition. The tracking 
of individual spots within the impact cloud shows that the westward jet at 56.5°S 

AW a , 

i 

latitude increases its eastward velocity with altitude above the tropopause by 5-10 ms' . 
The corresponding vertical wind shear is low. about 1 ms ! per scale height in 
agreement with previous thermal wind estimations. We found evidence for discrete 
localized meridional motions with speeds of 1-2 ms' 1 . Two numerical models are used 
to simulate the observed cloud dispersion. One is a pure advection of the aerosols by the 
winds and their shears. The other uses the EPIC code, a nonlinear calculation of the 
evolution of the potential vorticity field generated by a heat pulse that simulates the 
impact. Both models reproduce the observed global structure of the cloud and the 



dominant zonal dispersion of the aerosols, but not the details of the cloud morphology. 
The reflectivity of the impact cloud decreased exponentially with a characteristic 
timescale of 15 days; we can explain this behavior with a radiative transfer model of the 
cloud optical depth coupled to an advection model of the cloud dispersion by the wind 
shears. The expected sedimentation time in the stratosphere (altitude levels 5-100 mbar) 
for the small aerosol particles forming the cloud is 45-200 days, thus aerosols Were 
removed vertically over the long term following their zonal dispersion. No evidence of 
the cloud was detected 10 months after the impact. 
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1. INTRODUCTION 


On July 19, 2009 an object with an estimated size of ~ 500 m collided with 
Jupiter generating a debris cloud (thereafter called “the impact cloud” or IC) in the 
upper atmosphere with an east-west length of -4,800 km (Sanchez-Lavega et al., 
2010).. The debris cloud was subsequently dispersed by the Jovian winds at the impact 
latitude of 58.5°S (planetographic) (Hammel et al., 2010). The initial behavior of the IC 
globally resembled that observed for the individual SL9 impacnsyta- 1 994 (Hammel et 
al., 1995; West et al., 1995; Clarke et al., 1995; Sanchez-Lavega et al., 1995; Simon and 
Beebe, 1996; Banfield et al., 1996). However, we found some differences between the 
2009 and SL9 impacts. a 


% v ; ' 

For SL9, Sanchez-Lavega et al, (1998a) studied the long-term evolution of the aerosols 
caused by dispersion, mixing and vertical sedimentation. They found that after two 
years, a thin aerosol layer formed a wide band that fully covered the whole latitude 


circle, indicating dispersion in latitude by the action of the usually elusive meridional 


motions and low sedimentation rates consistent with small aerosol particles. The 
analysis of the dispersion of each individual element of the 1994 cloud debris (i.e. 

Ik 




produced by each impact) was complicated by the mixing of material between adjacent 

S\ 

impact sites. In contrast, for the 2009 event we had a single impact and the dispersion 
and spreading of the IC was easier to track, allowing for a better interpretation of the 
phenomenon. The study of the dispersion and evolution of the IC gives information on 
the physical state of the Jovian atmosphere allowing a direct inference of the vertical 
structure of the winds in the stratosphere, in particular at its upper parts where the 





aerosols can be used as tracers of the flow. It also gives information on the 
microphysical processes that took place within the cloud and led to sedimentation of the 
debris. Impact events such as this one and Shoemaker-Levy 9 in Jupiter allow us to 
compare the disperson and transport of energy and material left in the atmospheres of 
both planets, the induced local dynamics and aerosol microphysics (particle 


sedimentation, coagulation and coalescence). For the Earth, such studies are enabled by 
observations of cloud ashes produced during volcanic eruptions that inject aerosols in 
the upper troposphere and lower stratosphere (see e.g. Boville et jfc Mitff Stowe et al. 
1992; Trepte et al. 1993; McCormick and Veiga, 1992; Jonsson et al.. [996; Allen et al, 
1999). A comparison of these processes in the Earth ^itnfWpiter is enabled only by 
observations of these rare impact events. 

In this paper, we concentrate on the debris cloud evolution in 2009, from its very 
first detection on July 19 through December 31 when the cloud was barely visible. . The 


work represents a large effort by an international team that gathered images of the event 
over a broad wavelength range with a battery of telescopes and observatories, most of 
them listed in previous works on the impact (Sanchez-Lavega et al., 2010; Hammel et 
al., 2010; de Pater et al., 2010a; Fletcher et al., 2010; Orton et al., 201 1), as discussed in 

Section 2. We focus this study on images taken in two wavelength ranges: (1) visual 

W # # 

1 

range, broadly - 250 - 950 nm; and (2) near-infrared hydrogen and methane absorption 

jft V 

I bands (2.12 - 2.3 pm). Within each range we analyzed images taken with a variety of 


AL 


filters that have different bandwidths. The aerosol cloud can be well separated as a dark 
feature from the quiescent background atmosphere due to the high particle absorption in 
continuum wavelengths (Hammel et al., 2010; de Pater et al., 2010a), or as a bright 
feature on images taken in the gas absorption bands due to the high particle altitude. 


From the evolution and motions of the IC’s aerosol cloud elements reported in Section 
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3, we derive the global wind structure in the Jovian stratosphere at the impact latitude 
and discuss the presence of local regional dynamics. In Section 4 we present two 
dynamical models to explain the aerosol cloud evolution: a pure advection of the cloud 
by the winds and their shears, and non-linear EPIC code simulations that served to put 
an upper limit to the energy deposition in Jupiter’s stratosphere. Finally, from the 
reflectivity measurements and a radiative transfer model we calculate the aerosol optical 
depth evolution that is compared to that expected from the zonal dispersion and particle 
sedimentation. 

.-x-Jxv:- 

%# W 


2. OBSERVATIONS AND IMAGE ANALYSIS 



The images obtained in the visual range |f 250-950 nm) correspond to two main 
sets. The first set is comprised of images from July 19 to September 20 provided by a 
large team of contributors to the International Outer Planet Watch (IOPW) and the 
Planetary Virtual Observatory Laboratory PVOL data base (Hueso et ah, 2010). This 

data set has a daily temporal coverage of the IC evolution with small telescopes. The 

A ^ 

second set wgS the high-resolution observations obtained by the Hubble Space 

Telescope (HST) Wide Field Camera 3 (WFC3) during three runs in July 23. August 3 

A ' % 

and August 8, and September 23, previously analyzed by Hammel et al. (2010). 

Is., 

jm 


n % 

%:✓ 

The IOPW data set was obtained with a variety of telescopes, typically with 
apertures of 0.3-0. 6 m but with some observations obtained in larger telescopes with 
apertures of 0.81-m (IAC80, Canary Islands) and one set from the 1.0 m (Pic du Midi 
Observatory). At the small telescopes, fast imaging techniques (30-60 frames per 
second) with well-developed image processing techniques were used to improve the 
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resolution compared to single long-exposure (seconds) frame CCD cameras used in the 
larger telescopes (Cidadao. 2001; Grafton, 2003). Diffraction and seeing conditions 
limit the resolution to a maximum of about ~ 700 km/pixel for the best cases. Most 
images were sent to the IOPW-PVOL data base as processed color composites made 
with broadband filters in the range ~ 400-650 nm, and therefore they are not calibrated 
in intensity (Sanchez-Lavega et al., 2010). In addition to these “continuum wavelength 
composite images”, a small set of single images was obtained with a narrowband 
methane absorption filter at 890 nm. The IC is darker than surrounding clouds in the 
visible continuum, from the ultraviolet to the red region andsin thd color composites 
mentioned above, due to the intrinsic absorption of the4C "Slrdiols (Sanchez-Lavega et 
al., 2010; Hammel et ah, 2010; Orton et ah, 201 1 ; de Pater et ah, 2010a). In the 890-nm 
images, however, the IC is bright relative to the surroundings indicative of a relatively 
high altitude as previously indicated. These observations are important for two reasons. 


First, they provide an extended temporal coverage. Second, there are some instances in 
which images were made by different observers on the same day, allowing us to 
distinguish between small or diffuse cloud details and artifacts created by the image 
processing. The high-resolution HST imaging was used to complement the IOPW 

survey and to serve as a guide for the long-term cloud tracking of the IC features. They 

.# y '' 

were alsogused to measure the sizes of the small dark spots, as well as the edges of the 
zonally extended diffuse IC very accurately: the resolution was of the order of 200 
km/pixel. The details of these observations are presented in Hammel et al. (2010). 


The second large set we have analyzed for this paper covers the NIR spectral 
window from ~1 to 2.3 pm. The observations were performed by a large international 
collaboration involving larger telescopes around the world, as summarized in Table 1. 
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Further details for some of these observations were presented by Orton et al. (201 1) and 
de Pater et al. (2010a). We selected images taken in the NIR absorption bands at 2.12 by 
molecular hydrogen and at 2.23-2.3 pm by methane, both sensing the IC at higher 
altitudes due to the limited photon penetration in these bands. They allow us to 
distinguish the IC from the background even more than at 890 nm, because the IC is so 
bright relative to the surrounding background. 

All IOPW, HST and NIR images were processed and navigated ising the 
software LAIA based on the VICAR code (Cano, 1998). The methodoldfy employed to 
perform position measurements (longitude-latitude) of the features in time and extract 
motions from their tracking is detailed in Garda-Mclcndi&.et al: (2009). 

[Table 1 ] 


The NIR images were calibrated in intensity to determine the absolute 
reflectivity of the IC, as described in section 5. These reflectivity data were used to 
analyze the optical depth evolution of the IC during the observing period. Additionally 
we performed: The tracking of individual features in the images from July 20 to 

er this date, although large cloud condensations were still visible, 

'W 

their tracking was difficult and only the east and west “edges” of the extended IC 

-M, ' 


September . 



40 mm*-' 

> features were tracked. This was possible until October 29. Afterwards and until 
" December 31, the IC became homogeneous and largely extended and only area 
reflectivity measurements were possible. In summary, the reflectivity measurements 
cover the whole period from July 20 to December 31, 2009. The best coverage was 
during the first post-impact dates, from July 20 to August 13, when most observatories 
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contributed to the study. After this date, most images come from the IRTF-NASA team. 
Details on the exact dates of coverage are given in Table 1. 

3. EVOLUTION OF THE CLOUD DEBRIS AND WIND MEASUREMENTS 


The first image of the IC on July 19 (Sanchez-Lavega et al., 2010) and the 
higher-resolution HST images of July 23 (Hammel et ah, 2010) showed thPlinitial IC 
with two main components: a main dark streak surrounded in its westward flank by a 
halo or crescent that spans between planetographic latitudes -53. 5 & arid -60° (Figure 1). 
Planetographic latitudes (<p g ) are used along the text. Platietrigniphic and planetocentric 
( (p c ) latitudes are related by tan <p c = (R p / R E ) 2 tan® (Sanchez-Lavega. 201 0), being the 

polar and equatorial radii of Jupiter R P = 66854 km and Re = 71492 km (Lindal et al, 

1981). Afterwards, the IC expanded preferentially in the zonal direction by the 

4^. 

>x 

dominant zonal winds. Figure 1 shows a temporal map series, based on the IOPW and 

A‘ Wfr...#" 

gf ^ w 

HST images as observed in the visible spectral range, where the IC appears dark and 
increasing in size between July 19 and September 20. Figure 2 shows the temporal map 






series of the JC between July 20 and October 29 as observed in the near infrared 2.12- 


.,#■ .a % 




2.3 absor^jb^^ands (NIR, Near Infrared) where the IC appears bright. The same 
genital rnbrphology of the IC is observed in the images taken in the 890-nm methane 
absorption band. A simple comparison immediately shows that the same features are 
observed within the IC in the continuum and methane bands but with reversed 
brightness. The identical cloud patterns within the IC observed at these wavelengths 
suggest either (a) both visible and near-IR wavelengths sense the same altitudes; or (b) 
the debris material is distributed uniformly over a range of altitudes (de Pater et a. 
2010a). In this paper we assume that we detect the IC at the same altitude levels in the 





visible and near-IR. Since the debris is detected in the 2.3|im absorption bands, the IC 
must be located at and above the tropopause (Hammel et al., 2010; de Pater et al., 
2010a). Thus the differences in darkness (absorption in the visible continuum) or 
brightness (reflectivity within the methane absorption bands) within the IC mostly 
indicate differences in optical thickness, but not large altitude differences. 


[Figure 1] 


[Figure 2] 






% 


Just as in the observations of the SL9 impact in 1994, the clouds spread 
preferentially in the zonal direction i.e. that of the prevailing winds, indicating that the 
aerosol transport by the zonal winds dominates the evolution of the IC morphology. 
However, the cloud is not homogeneous in its reflectivity (i.e. in its optical thickness) 
but shows a distinct distribution of clumps or spots with higher aerosol density, 
particularly evident during me first month after the impact. These spots were used as 


tracers of the motions, measuring their longitude and latitude versus time. But 

Jr \ 

importantly, as the debris cloud was dispersed by the winds, it also showed an intrinsic 
evolutiomwith new spot formation and in other cases with some of the spots dividing 




. into two or more elements (Figures 1 and 2). This aspect wi 


will be discussed in section 


3.3. 


3.1 Longitude tracking 




The evolution of the IC was determined by direct latitude and longitude 
measurements of the features on the navigated images and on cylindrical map 
projections. The tracking of the condensation spots was performed by visual 
identification of the same feature on successive images. We separate the tracking into 
three main sets, two corresponding to the zonal drift of the features and one to t&eir 
meridional motion. To determine the zonal drift rate, we used two methods. In the first 
method, we measured the East and West longitude limits (the "edges").©# thellC in the 
two wavelength ranges (visible and NIR) as a function of time. Tliis is shown in Figure 
3. The IOPW-HST tracking data cover the period July 19 - Aug. 14 and the NIR 
tracking data cover the period July 20 - Oct. 29, 20^.^ Afterwards, only brightness 
studies were performed until December 31, 2009, after that time the planet was difficult 

vX; . 'A 

to observe due to its proximity to the Sun (solar conjunction took place on February 
28.). The expansion rate of the IC edges can be reasonably fitted by a linear function 
L{III) = 0 )At + with l Jill) the System III longitude, co the zonal angular velocity 

is the reduced Julian Date, where J.D. = 2455032 


x .# 

.•vs 


(deg/day) and At = J.D. - 2 .. „ 
for July 19 (at 12 U.T.): 




Eastern (“leading”) edge: L(IH) = 0.51* At -2259 


Jr* 


# 


Western ("trailing") edge: I Jill) = — 1 .02 • At + 543 1 




II 

; :V 

1 


( 1 ) 


( 2 ) 




According to these edge displacements, the prediction is that the east and west edges of 
the IC, which moved in opposite directions, should have encountered one another at 
longitude ~ 67° and on A/ ~ 5258 that corresponds to early March, i.e. during the solar 
conjunction. After solar conjunction, similar NIR observations were made in May, 



2009, but by then the reflectivity at the IC latitude was indistinguishable from that of 
the pre-impact atmosphere. 

The averaged zonal velocity, <u>, as retrieved from these linear fits have low 
errors a(< u>)~ Ax/ At that for positioning precision Ax <1000 km and tracking time 

At ~ 60-100 days give cr(< u>)~ 0.5-1 ms ! . 


[Figure 3] 

In the second method, we tracked the longitude-latitude positions of individual 
spots within IC in the visible continuum (July 19 - Sept. 20) and in the near infrared 
(July 20 - October 29). The longitude drift charts of the individual tracked features are 
shown in Figure 4 for the visible and in Figure 5 for the near infrared. Each tracer was 
typically measured at least in 5 different images at different times. The minimum time 
interval to follow a tracer is 3 days and we found that a single tracer does not usually 
survive more than about 15 days, or else can not be identified as the same feature. In 
total, we tracked 53 spots. For completeness, we added to the tracking chart those points 
previously reported from our HST measurements (Hammel et al., 2010). It must be 

noted that -some tracers were measured simultaneously in the IOPW, HST and NIR, but 

#’ 

usually spanning different time intervals, making the results consistent. The longitude 
drift of each feature was fitted to a linear function obtaining the zonal angular velocity 
of the feature in System III and the zonal velocity («,). As before we expect from the 
tracking zonal velocity errors for each individual feature cr(w,) < 1-3.5 ms 4 
(positioning precision Ax < 1000 km and typical tracking times At > 3 - 15 days ). 




[Figure 4] 


[Figure 5] 


3.2 Latitude tracking 


% 


Although the spreading of the IC occurred mainly in the zonal direction due to 
the dominant zonal winds and its meridional shear, a meridional: expansion for some 
features was also detected. Measurements of the HST images on August 3 and 8 showed 
that whereas the northern limit of the IC remained at the initial latitude of -53.5°±0.5°, 
the southern limit was at latitude -61.4°±0.5°, i.e. closer to the pole than initially 
observed. This indicates that a poleward drift of the IC occurred during the first 2-3 
weeks at an approximate rate of 0.1 "/day. Our long-term measurements of the 
meridional limits of the IC (July 19 - December 30) indicates that the IC aerosols 
remained confined between these two latitudes -53° and -61° where two eastward jets 
have their peak speeds and where the gradient of the potential vorticity becomes large 


(see sections 3.3 and 4.3). 




however. we detected a faster meridional motion, both in the visible and NIR 


images, for a single spot that formed in the first days after the impact (Figure 6). 
^Tracking its motion (longitude and latitude) from July 23 to August 5 revealed that it 
translated from latitude -56.0° to -53.5°. Navigation errors in the high resolution HST 
and NIR images translate in latitude uncertainty measurements of 0.2°. The total 
meridional displacement of the feature was ~ 2.5° (equivalent to 3120 km) between July 
23 to 30, equivalent approximately to a meridional drift rate of +0.3°/day. In summary, 




mvmmmv 
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IC aerosols were observed to drift both poleward and equatorward, specifically during 
the early stages of the IC evolution. 


[Figure 6] 


3.3 Wind measurements 


The measured longitudinal drift of each individual clou<f 'element was 
transformed to the zonal wind speed at its mean latitude. Iri Figure 7 we present the 
zonal wind profile that we derived for all the features %e flteasured within the IC and 
compare it with one measured using regular Jovian features using HST images obtained 
two years before the impact in 2007 and intbe same HST images in 2009 used to track 
the IC features but this time using discrete features outside it. A parabolic fit to the IC 
velocity points is shown in Figure 7 (right) to serve as comparison with the background 
HST profile and to have an idea of the velocity dispersion of the IC elements. Two 

results emerge from these profiles. First, it is evident that the wind velocities of the IC 

'm.. 

elements are h|| than those of the reference background wind profile. Since the 


aerosol 



is, placed at a higher altitude than the background cloud (Hammel et al., 


20 If; de bater et al., 2010a), these measurements indicate that the winds increase with 


w ,,„. 


| altitude within the latitude band where the impact aerosols formed (latitudes -53° to - 

m _ -yy 
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61°). The speed increase is within 5 to 10 ms' 1 (which is above our estimated 
measurement errors). According to Hammel et al. (2010). de Pater et al. (2010a) and 


Orton et al. (2011), the cloud debris extended vertically from 1-20 mbar to 200-300 
mbar (about 4-5 scale heights H; 1 H = 22 km). Assuming that the background wind 
profile is determined from cloud tracking at the 500-mbar pressure level, the vertical 
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shear of the zonal wind is expected to be of the order of +1 ms -1 per scale height which 
is in agreement with the thermal wind speeds derived from Voyager- IRIS and Cassini - 
CIRS measurements at these latitudes (Simon-Miller et ah, 2006). 

[Figure 7] 


Second, all the tracking points show a significant velocity dispersion at a given 
latitude with values varying by up to 10-15 ms 1 (e.g. at latitude -56°, Figure 7). This 
dispersion is well above the time-averaged velocity measurement error, which we 
estimated before to be <r(w,.)< 1-3.5 ms' 1 from the large temporal tracking of each 
feature. Since we estimate the averaged error in the IStltudinal position of each feature 
to be ± 0.5°, this indicates that although the features followed the mean atmospheric 
flow as nearly passive tracers of the zonal winds, they also showed intrinsic local 
motions. In particular, there was evidence of clump formation (condensation spots) 

a $ 
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within the IC during its expansion with clump size of order 500 - 1000 km (Figure 8). 


Similar clump formations and fast initial motions were observed in some of the ICs 




generated by 


impacts (Sanchez-Lavega et al., 1995; Simon et al., 1996). The 




possible nature of this clump formation and their motions within the IC will be 

jm 

disctissed below. We also note that this velocity dispersion is similar in magnitude to 
the fluctuations in zonal speeds as measured on timescales ranging from a Jovian 
- rotation up to weeks or months (Asay-Davis et al. 201 1). 


Finally the latitudinal drifts of some features that were measured during the first 
two weeks after the impact were converted to meridional wind velocities. At latitude - 
60° the motion was poleward with a velocity of 0.8 ms’ 1 whereas in the middle of the IC 
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at -55° the motion was equatorward with a velocity of 2 ms' 1 . As indicated above, no 
further meridional spread was seen, and the IC remained confined between latitudes - 


53.5°±0.5° and -61.4°±0.5° until the last observation in which it was detectable. 


[Figure 8] 


3.4 Reflectivity measurements in the NIR absorption band 



In order to study the aerosol content evolution withitl the IC. we performed a 
long-term measurement of its brightness reflectivity I/Fin the NIR methane absorption 
band at 2.2 - 2.3 pm. We used for this study the images coming from the NIR dataset 
which corresponds to the contributing observatories numbered (3) to (9) in the first 
column of Table 1. According to methane bhnd coefficients derived by Karkoschka and 
Tomasko (2010), the optical depth for the gas absorption at these wavelengths is unity 
at pressure levels ~ 1-10 mbar; thus, we are sensing the top of the IC. The reflectivity 
was measured relative to that of close latitudes, not affected by the impact, to account 


for the geometric variability across the Jovian disk, most importantly the center to limb 
variation. An absolute reflectivity calibration was performed by reference to the 

JP TV 

reflectivity curves in de Pater et al. (2010a). Since the cloud brightness structure was 
not spatially homogeneous, in particular during the first two months (Figure 2), an 
veraged IC reflectivity was obtained for each date from measurements across different 



points of the IC between July 19 and September 20. For most of the data points, the 
geometrical configuration corresponds to the cosine of the incident and emission angles 
with values n ~ juo ~ 0.6 (corresponding to the central meridian at about -58° latitude). In 



Figure 9 we show the measured data point reflectivity for the IC as a function of time in 
this methane band together with a simple decaying exponential fit 


-(0 = 0.0025 + 


0.017 exp 


f ( At - 5040) ' 


(3) 


where At is the reduced Julian date (Julian Date - 2,450,000; Julian Date of 19 July 
2009 12:00 UT = 2,455,032) and to is the e-folding timescale for the reflectivity decay 
which we found to be ~ 15 days from the fit. 




[Figure 
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4. ANALYSIS OF THE DEBRIS CLOUD EVOLUTION 
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In order to interpret the structure and evolution of the IC we must first review 

„ - 0 ' 

the known dynamics otfhis Mpiter region. As shown in Figure 7, the impact occurred in 


a band where the flow is weakly westward, bounded on the north and south by two 

# % 
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eastward jets With peak speed of - 30 ms at latitude -61° and ~ 45 ms at latitude -52° 

(Garcia-Melendo and Sanchez-Lavega, 2001). The southern polar region is 
characterized by the presence of planetary scale waves that manifest in the morphology 
of the high altitude hazes and cloud field as a wavy pattern that encloses the polar 
region (Sanchez-Lavega et al., 1998b; Vincent et al., 2000: Barrado-Izagirre et al„ 
2009). There are different waves whose latitudinal location depends on the observing 
wavelength. The contrast is largest in the ultraviolet where they appear dark due to haze 
absorption. They are bright in the methane absorption bands (890 nm, 2.23 - 2.3 pm) 
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where they vary in brightness, following a trend similar to that of the IC. Sanchez- 
Lavega et al. (1998b) and Barrado-Izagirre et al. (2009) demonstrated that the observed 
dynamics is consistent with Rossby wave activity, based on their motion and dispersion 
relationship. These waves serve as a diagnostic tool to test the atmospheric response to 
the impact at nearby latitudes. -<$^4 


Figure 10 is a series of map projections of the southern polar regi6h ; :at these 
wavelengths showing the IC and the south polar waves. The impact occurred midway in 
latitude between two waves: the poleward main wave that is prominent in the methane 
band images centered at latitude -65°, and an equatorward wave that is similarly 
prominent in the ultraviolet and is centered at -48°. These waves are prominent in the 


HST UV images at wavelengths 225 - 378 nm (Hammel et al., 2010) and subtly 
detected in thermal- IK imaging of stratospheric CH4 emission near 7.9 pm by VLT 


(Fletcher et al., 2010). 


ring the amplitudes and structure of these two waves in 


the HST images obtained in July 23. August 3 and 8. it seems that they were not 
affected at all by the : impact. The wavy pattern was always present at the impact 

longitudes, without apparent significant changes in the meridional amplitude and zonal 

0 / . ■■■0 

wavelength, indicating that the impact did not affect the dynamics of these permanent 
waves. The response of the atmosphere to the impact was restricted in latitude to the 
impacted band that had a meridional width of ~ 8°. 


[Figure 10] 


4.1 Global dispersion of the cloud debris by the zonal winds 




The observed dominant zonal spread of the IC and the exponential decay of its 
2. 23-2. 3pm reflectivity were originated by the aerosol dispersion and subsequent 
decrease in the particle number density and optical depth as a result of various 
processes. The most obvious is the dispersion produced by the prevailing zonal winds 
and its shear (meridional and vertical). The typical temporal scale for the aerosol 


dispersion due to the zonal wind shears is t n ~ (du/dy) ' ,(du/dz) ’ - 1-5 _ days, 

according to the measurements reported in section 3.3 and Figure 7. This ig a factor 3- 
15 faster than that indicated by the e-folding value for the measured reflectivity decay, 
suggesting that dispersion by wind shears was the main operating mechanism. 

To test this, we have run a simple adveetion model in which the IC is dispersed 
by the zonal winds and their shears. The model solves the adveetion equation for a 
passive tracer under a purelly zonal wind which depends on latitude and height. In this 
simple approximation no diffusion is incorporated and the only equation to solve is: 


dq dq n 
at ox 
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where q is the density of thO debris material and u is the zonal wind depending on the 
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latitude and hcight.The numeric code is a straightforward simplification of the two- 
dimensionat code used by Hueso et al. (2002). In this case the model domain consists of 




a 


• boraputed 




^ (30 grid points for longitude and latitude and 11 vertically detached layers are 

simultaneously to incorporate effects of a vertical wind shear. The horizontal 
grid points are separated by 40 km while the vertical layers are representative of close- 
in levels in the stratosphere. For the zonal wind we impose the parabolic fit to the zonal 
wind profile shown in Figure 7 in the central layer, with the winds vertically varying 
from one layer to the next by 1 ms" 1 faster or slower when above or below the central 
layer. Therefore the model incorporates a total vertical wind shear of 10 ms" 1 with each 



vertical layer decoupled from the rest. We started the simulations by placing a double 
oval-shaped cloud that mimics the main streak and the halo or crescent area of the 
observed IC (Sanchez-Lavega et al., 2010), oriented in the longitude-latitude plane as 


observed (Figure 11) and with the same structure in all 11 vertical layers. The cloud 
density is uniform within the two areas that compose it but differing in their value from 


the streak to the halo in a grey scale that represents the observed initial rafl|ptrl|ty. 
Figure 1 1 shows snapshots of the simulations with the cloud advected apd dispersed by 
the winds for a time period of 40 days. The simulation produce ahPi^vPlving debris 
cloud whose vertically integrated density decreases with an e-folding time of 10 days 
produced by the combination of both zonal and vertical Sheafs. Stripes in the figure 
appear from the limited number of vertical layers used which each vertical layer moving 
with a slightly different zonal wind compared to the others. 

Although the zonal expansion and growth of the dense parts resemble those 


observed for the IC (Figures 1 and 2), there are still differences with what is observed, 
as for example the faiiure (o reproduce the north-south “horn-shaped" form of the 
simulated expansion, that is hot observed. 


[Figure 11] 




To explain the differences between the simulated and the observed evolution of 
the global IC expansion, we must invoke other mechanisms than simply the structure of 
the zonal winds. Among them are the action of local thermal winds that could have 
developed from the temperature anomaly originated by the energetic impulse of the 
impact. From the wind data at this latitude the geostrophic balance can be assumed 
since the Rossby number Ro-u! f L~ 0.005 (here u = 10 ms' 1 , L ~ 7000 km. 




approximately the westward jet width, and / = 2Q . sin (p = 2.88xl0" 4 s 4 ). Orton et al. 
(2011) measured temperature differences of ~ 3 - 6 K between the impacted area and 
surroundings a few days after the impact. Assuming that this temperature contrast 
occurred over the same pressure level, the induced local thermal wind should reach a 
speed difference between pressure levels Po and Pi of the order of (Sanchez-Lavega, 
2010) A> ' ' 





the surrounding IC aerosols. The spots were concentrated in the anticyclonic part of the 
wind profile and their total number was about 14 on August 3. They disappeared in a 
matter of days probably by mergers between them, whereas others were dispersed by 
the wind shears, mixing with the background aerosol field. 


Their distribution within the IC and their structure suggest that these spots may 
be vortices resulting from local instabilities. The high resolution HST images |Hammel 
et al., 2010) do not show well defined spot boundaries, but instsad iiey seem to be 
elongated with irregular edges and filaments. The ambient fldw: in:the anticyclonic part 
of the IC has a local vorticity £ = du/dy = -4.2xl0' 6 s' 1 , about two orders of magnitude 
lower than the planetary vorticity / = 242 sin <p = 2.9x10 1 s' 1 . The interaction of the 
background flow with the spots will produce a small tangential speed at its periphery of 
the order V T ~ gR 12 = 1-2 ms" 1 , with R its size (effective radius). However, these spots 
could have other vorticity sources and thus be highly speeding rotating vortices. 
Unfortunately the HST resolution does not allow us to track features rotating around 


them. 
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In order to constrain the action of possible instabilities in generating these 
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clumps we can compare their size with some characteristic dynamical scales that are 

, # g * 

related to different flow instability mechanisms. For example, one important scale for 


the formation of baroclinic eddies in quasi-geostrophic (QG) theory of a continuously 
stratified atmosphere is the first (internal) Rossby deformation radius L D = NH/f (see 


m 

e.g. Gill, 1982; Vallis, 2006) where a factor 7t in the denominator has been omitted, as 
usually done. This scale has been shown to be important in models of regular Jovian 
vortices (Shetty et al., 2007; de Pater et al., 2010b and references therein). Here N is the 
Brunt-Vaisala frequency that in the Jovian stratosphere is - 1.9xl0‘ 2 s" 1 and H is usually 
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taken to be the scale height ~ 22 km, which gives L D ~ 1400 km. Simple scaling will 
give baroclinic eddy velocities (Vallis, 2006) v e ~{LJL d )<u>~ 4-8 ms" 1 where we 
have taken an eddy scale L e ~ 500 - 1000 km (that of the clumps) and a mean flow 
velocity <u> ~ 10 ms" 1 . The corresponding eddy time scale is T e ~(L e /v e )~ 1 day. 
Weaker vortices are obtained if, instead of this value for H, we use the vertical extent of 
the IC, i.e. we assume that the spots extend vertically along the layer occupied, by the 
IC, then Hie = H In (300 mbar/1-10 mbar) - 3. 5-5. 5 H, L D ~ 5000 - 8000 km and v e ~ 0.5 
-1ms . 

Since the vertical shear of the zonal flow is low, another important scale to be 


considered is that associated to the barotropic instability of the westward jet where the 
IC formed. Ingersoll et al. (1981) showed that Jovian westward jets violate the Rayleigh 
- Kuo barotropic stability criterion [5- d 2 u/dy 2 >0 being/f = df l dy = 2Q.cos(pl R p 
the gradient of the Coriolis parameter / (Sanchez-Lavega. 2010). The characteristic 
horizontal scale associated Jb tlis : i.ltistability is ~ itL p = xju/ j3 with growing times 

t b ~ 10 1 Juft (Kuo 1 949. 1973 , see also Stone, 1976). For u ~ 10 ms" 1 and B- 2.9x10" 

% 

12 ms" 1 at this latitude we get Lp = 1800 km and Tp ~ 20 days. The size of the clumps are 

A % 




within the order of magnitude of L D and Lp . This simple scale analysis indicates that 

jf - ' M 

baroclinic or barotropic instabilities could plausibly have played a role in the clump 

formation. 

' 
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4.3 Numerical simulations of thermally induced instability in the stratosphere 


We have studied the dynamical response of the Jovian atmosphere with the 
properties of the impact latitude to an instantaneous and localized thermal pulse using 





the EPIC code (“Explicit Planetary Isentropic-coordinate”, Dowling et al., 1998). Our 
purpose is to see how the potential vorticity (PV) field resulting from an impact and 
assumed to represent a passive tracer of the flow, mimics the shape, extent, structure 
and evolution of the IC. Details of EPIC model simulations of the Jovian atmosphere 
have been extensively described in our previous papers (e,g. Morales-Juberfas et al., 
2003; Garcfa-Melendo et al., 2005; Legarreta and Sanchez-Lavega, 2008; Sdnchez- 
Lavega et al., 2008). 




We introduce a localized thermal pulse in an EPIC model atriiosphere that tries 
to simulate the thermal anomaly produced by the impact. We have run series of 
numerical simulations to draw maps of the evolution of the potential vorticity (PV), 
selecting two atmospheric levels for the representations, located at ~ 45 and 400 mbar, 
that corresponds to the IC formation altitudes. 


The model atmosphere depends on several free parameters, the most important 
being the vertical and; meridional structure of the zonal wind velocity U(y,z) and the 
vertical thermal structure represented here by the static stability or Brunt-Vaisala 


frequency. The three-dimensional structure of the wind field is separated as a product of 

jm 

two functions U(y, p) = u h (y)u v (P) where u h (y) is the wind profile at the cloud tops 






y; and u,(p) is a nondimensional vertical amplitude factor. For the meridional wind profile 




u h (y) at cloud-top level we used that derived from HST images in 2007 (Figure 7) 
locating it at a pressure level Po = 500 mbar. For the vertical structure we used a two- 
segment linear dependence in InP as in Garcfa-Melendo et al. (2007). For P 0 < 500 mbar 
we tested winds that slightly decrease, are constant or slightly increase with altitude 


according to the expression u v (P) = 1 + m In (P/ P 0 ) with an e-folding scale m in 
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steps from +7, 0, -10 scale heights, i.e. winds that change their speed in 1/e in these 
scale height values. This range covers the observed vertical wind shear at the impact 
latitude (i. e. a westward jet that increases its eastward speed with altitude) as well as 
other possibilities to test the sensitivity of the atmospheric response to the vertical shear 
sign. For P > 500 mbar, we assume that the winds remain constant with depth. For the 
thermal structure, for pressures lower than Pq we used the thermal profile obtained from 
radio occultation experiments (Lindal et al., 1981); below cloud tops (pressuTfes greater 
than Po), the thermal profile was extrapolated by a wet adiabat characterized by the 
value of the static stability of the lower levels (see details in Garci'a-Melcndo et al., 
2005; Legarreta and Sane hez-La vega, 2008; Sanchez-Lavega:ef al.. 2008). 


Two groups of simulations were run m a channel that spans latitudes from -38° 
to -68° using 8 vertical layers: (1) A low-resolution set with a model domain that covers 
60° in longitude and with 256 x 128 grid elements in longitude and latitude respectively 
and (2) A high-resolution sltjvitHl model domain that covers 60° in longitude and 15° 


„ -m m 


in latitude with 512 x 128 pid elements In this second set, the spatial resolution was 




0.12° per grid element equivalent to 149 km in latitude and 80 km in longitude (the 


range in lphgitude for the domain we used is 60 km to 100 km). We performed a total of 
160 simulations each of which extended between 25 and 50 days in time. The 

disturbance produced by the impact was simulated by introducing an instantaneous 

% ,# 

initial heat source in the EPIC code with a power ranging from 10 to 10,000 Wm‘, and 

%f "" 

occupying a size of 10° in longitude by 3° in latitude (approximately that of the initial 
IC, also similar to the temperature anomaly in Orton et al., 2011) and extending 
vertically from the upper troposphere to the lower stratosphere (pressure range from 1 


mbar to 1 bar). We used an hyperviscotity coefficient (Dowling et al., 1998) v 




0.57xl0 27 mV 1 that for a scale AL - 400 km gives a diffusion time scale 
A t D = (AL) 6 Iv 6 - 100 Earth days which is above the length of our simulations. 


The results, in the form of the evolution of potential vorticity (PV) maps at the 
two selected levels (specify levels), show that the most important parameter controlling 
the PV-field is the impact power. For input powers ~ 500 Wm~ 2 the PV-field evolution 
follows a pattern that resembles one produced by the passive advection of the IC as 
discussed in section 4.1. The simulated PV field reproduces the essentials of the IC 
evolution, i.e. its zonal expansion along the anticyclonip, side of the zonal flow. In 


Figure 12 we show the result of one of the best simulations for which the PV field 
evolution at the 1 80-mbar altitude level most closely resembles the observations of the 
aerosol reflectivity evolution. When higher input powers are introduced (above 500 
Wm' 2 ) the PV-field tendency is to form a vortex rotating in anticyclone sense following 


the ambient wind shear vorticity. The vortex tends to produce wave instabilities and 
turbulence around it and in the adjacent latitudes (those corresponding to the polar 
waves discussed in tie introduction of section 4), that are not evident in the 
observations of the cloud and thermal fields. The simulations produce spots of PV 


within the large vortex, rotating around its centre. One possibility for the clump 

£ i 

foril^tipri is that they result from the interaction of the IC with pre-existing, underneath 
to the IC, regular Jovian meteorological features (e. g. vortices and waves). Our EPIC 
model cannot capture this since these meteorological phenomena were not introduced in 
the initial simulations. 


The potential vorticity field as drawn from the EPIC calculations (Figure 12) 
shows large values of its gradient at latitudes -52° and -61°. This potential vorticity 
gradient could have acted as a barrier to the meridional transport of the IC aerosols (see 



section 3.3), as seen for example in the compound transport in the Earth's stratospheric 
ozone hole (McIntyre, 1989). 


[Figure 12] 


4.4 Aerosol optical depth evolution and sedimentation. 


In order to reproduce the observed variation of reflectivity in the 212-2.3 pm 

methane absorption band, we used the same vertical structure ofjhe impact debris as 

was used in previous studies for dates closer to the impact (de Pater et al., 2010a). In 

that work, the impact site reflectivity was fitted from thfe ; ,ulHs^iblet to the near infrared 

with essentially two layers of aerosols. The upper layer contains particles with a mode 

radius a = 0.75 pm, a particle density of 1.2 part/env\ and with a vertical range of 10 - 

1 10 mbar; the lower layer has particles with the same mode radius, but a particle density 

of 5.6 part/cm 3 and distributed in a higher pressure range of 1 10 - 200 mbar. Since the 

contribution of the layers below 200 mbar to the observed reflectivity at 2.3 pm is 

■ 

negligible (see contribution functions at Fig. 14b from de Pater et al., 2010), we fitted 


the observed 







ral dependence of the reflectivity by varying the particle optical 




thickness at theonly two aerosol layers above this altitude in de Pater et al. (2010) 
model. The only free parameter in the model was therefore the particle concentration at 
such layers. We assumed three different scenarios for the density evolution: (A) 
A.-::;.. Particles from both layers were removed simultaneously and at the same rate; (B) the 

1§f 

upper particle layer was emptied first, followed by the lower second layer; and (C) the 
lower layer was emptied first, followed by the upper one. Dynamically, scenario (A) 
seems more realistic, but scenarios (B) and (C) provide upper and lower limits for the 
particle concentration retrieval at each layer. 
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Results are shown in Fig. 13. Both layers have to be drained to explain the 

observed reflectivity decay (Figure 9), and all three scenarios provide a similar decay 

law. The second lower layer (which was also initially denser) is about one order 

magnitude denser at the final stages, with a particle concentration lower than ni ower <0.5 

, 

part/cm , whereas the first upper layer reaches densities not higher than 0.03 

vxx ; ;., 

part/cm\ Scenario (B) provides a fast depletion of the upper layer and a much slower 
depletion of the lower one. just the opposite of what scenario (C) 


m. 


A. 


[Figure 13] 


According to this model, the total particle optical thickness decay within the IC 
can be fitted by an exponential function of the form 


X.v:;.. 


T(t) = 0.014 + 


0.41 ex 


, % J 

t D 


(5) 


where we find t a = 14.5 days, similar to the 15 days that we found from our fit to the 
observations (eq. 3). This decay law is also similar to that observed in the long-term 
analysis of the SL9 impact clouds whose evolution was more complex due to the 
•^^^Jtnergers between close but different impact clouds and their subsequent dispersion 
(Sanchez-Lavega et al., 1998a). For the particle concentration at each layer we get 


n upper (0 = 0-028 + 


0.84 exp 


f (At -5046)^ 


( 6 ) 





n iower( 0=0.13 + 


4 exp 


( (Af-5046)^ 


l D J 


(7) 


where time is given in days in the reduced Jovian date scale (section 3.1) and the 
particle concentration is given in part/cm'. The exponential decay law is consistent with 
a horizontal spread of the aerosols by the wind shears as described in Sanchez-Lavega et 
al. (1998a). Assuming that the total number of particles in a giveq .yoluntehA remained 

constant during their transport, the column density n and the area of the site A obey the 

W 

continuity equation 


1 dn 
n dt 


}_dA 
A dt t n 


( 8 ) 
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whose solution is equation (6)-<7) A similar relationship can be formulated for the 
optical thickness t = with Q the scattering efficiency factor, with a solution 




of the type of eqaiilion (5). 


Particle sedimentation can also contribute to the evolution of the IC optical 

£ f 

f h. We have studied the particle density evolution due solely to sedimentation as in 

Its- M 

Sanchez-Lavega et al. (1998a). The aerosol sedimentation time depends on the particle 

’ll# . 

size and atmospheric properties that can be characterized by two dimensionless 
numbers, the Knudsen number (Kn) and the Reynolds number (Re). In the Jovian 
stratosphere (altitude levels 10-20 mbar) we ha \eKn = £/a~ 2 and Re = (2 paw) It ] ~ 




1.3x1 O' 6 , with £ the gas molecular mean free path, a - 0.75 pm the particle radius, p the 
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gas density, w the vertical fall velocity and jj the dynamic viscosity. The values for the 
latter parameters have been taken from Sanchez-Lavega et al. (1998a). These values 
bound the dynamical mechanisms that act on the aerosols between those of the laminar 
and gas kinetic regimes, giving a sedimentation time for the particles t s ~ 45 - 200 days. 
This suggests that the aerosol sedimentation is only important over long time scales, i.e„ 
removing aerosols via sedimentation well after the IC debris has been dispersed by the 
winds, thus effectively after November 2009. 


5. CONCLUSIONS 




The evolution of the IC’s debris from the Jupiter impact on 19 July 2009 mimics that 


seen for the SL9 impacts in 1994, although in 2009 it was easier to determine the 

evolution as it was generated by a single impact. The structure and long-term evolution 

m, |? 
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of the IC showed that in the stratosphere (up to ~ 1-10 mbar altitude level) and within 


the latitude band 
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by the IC (from 53.5°S to 6S1.5°), the zonal motions were 


still dominant, although they were decreasing their eastward motion. 
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I The westward flow at the impact latitude became weaker with altitude, i. e . , it 

increased its eastward velocity upward above the tropopause. The two bounding north 

iL M 

and south eastward jets appeared to be less intense in their velocity peak. Thus globally, 


the zonal flow became weaker with altitude. This confirms previous studies on the wind 
vertical structure based on the temperature measurements in the stratosphere and 
application of the thermal wind relationship (Simon-Miller et al., 2006). 





In addition, this study confirms that the meridional motions in the stratosphere 
were also weak and, at least at this latitude, confined between the two contiguous 
directionally opposed zonal jets. This is consistent with the thermal-IR analysis (Orton 
et al., 2011; Fletcher et al., 2010), which also indicates horizontal divergence over the 
impact latitude. However we note that this is at odds with the SL9 impact debris 
meridional motions that spread much more in latitude, moving northwards, j®j5howed 
by Sanchez-Lavega et al. (1998a) and by HCN measurements (Griffith el al., 2004 ). 


The evolution of the IC was mainly produced by the%erosol dispersion by the 
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zonal wind flow, with the particles advected according,td)£he wind shear. However, 
explaining the details of the IC structure, such as the formation of clumps (condensation 


spot), is not easy and deserves further study. These clumps may be vortices generated 
by the wind shear under the action of some kind of dynamical instability. Local motions 
induced by the temperature anomaly produced by the impact, superimposed to those 


related to the Jovian background dynamics, could have also contributed to the IC 


evolution. Because of the Jovian stratospheric properties and small particle size 

^ % 

composing the IC ’•the residence time of the aerosols is large, with sedimentation acting 

J? a . 

effectively in removing them after several months. This is in agreement with the 

■yisf If 

J0- 

observations that showed the aerosol particles still present almost 6 months after impact 

Ik .# 

and undetectable 6 months after the impact. 

IL JF 




* : 


Since we did not detect changes in the system of circumpolar waves located 
north and south of the impact region we conclude that the dynamical effect of the 
impact remained confined in latitude. This indicates that the response of the atmosphere 





to the energy deposited by the impact was merely local, not affecting the dynamics of 
the closest latitudes. 
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TABLE 1. Observatories contributing with images to this study 


Observatory 

Telescope 
( diameter ) & 
Instrument (*) 

Wavelengths 

Date coverage 

Month (days) % 

IOPW 
network (1) 

(1) 

350 - 950 nm 

Jul (21,22,23,24,25,26,27,29,30,31:), 

Aug( 1 ,2, 3, 4, 5 ,6, 7, 8, 9, 1 0, 1 1 .12.13,14,1 5,17, 
19,20, 22,24,26,27,29,30), 

Sep (1,5,6,8,10,1 1,13,15,17,18,20) 

HST (2) 

WFC3 

234 - 924 nm 

Jul (23), Aug (3, 8), Sept (2H ; 23) 

Calar Alto (3) 

3.5 m - 
Omega 2000 

1.6-2.17 pm 

Jul (20,22) 

ESO-VLT (4) 

8 m - NACO 

1.6-2.48pm 

Jul(26) 

TNG (5) 

3.6 m - AdOpt 

l-2.3pm 

Jul(27), Au||3.4;!3) 

WHT (6) 

4.2 m - Ingrid 

l-2.3pm 

Aug (4,% : P 

Keck (7) 

10 m- NIRC2 

1.6-2.12 pm 

Jul (20fe24:25). Sep. (7,10) 

IRTF-NASA 

(8) 

3m- SpeX 

2.12-2.23 pm 

jgjj % 

Jul (2024.26) 

Aug (2,4,13,14,15,16,18,19,23,25,30) 
Sep (2,7,8,9,14,20,23,25,30) 
lOct (1,19,20,25,27,29) 

Nov (44,8,18,22,23,24,25) 

Dec (4,5,29,30,31) 

May (14, 28) [2010] 

TCS (9) 

1.5 m - CAIN3| 

l- 2.3pm 

Oct (7, 23) 


M 


Notes: 

(*) Telescope diameter and instrument/detector employed. Details for the 
instrumentation can be obtained in the websites of the Observatories. Other details 
given in de Pater et al. (2010) and Orton et al. (2011). 

( 1 ) fOPW (http://www.ehu.es/iopw) and database PVOL (http://www.pvol.ehu.es/) . See 
Hueso et al. (2010). 

(2) Hubble Space Telescope. See details of the observations in Hammel et al. (2010). 

(3) Calar Alto Observatory, Centro Astronomico-Hispano Aleman (Almerfa, Spain) 

(4) European Southern Observatory - Very Large Telescope (Chile) 

(5) Telescopio Nazionale Galileo (Canary I., Spain) 
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(6) William Herschel Telescope - ING (Canary I., Spain) 

(7) W. M. Keck Observatory (Hawaii, USA) 

(8) Infrared Telescope Facility-NASA (Hawaii, USA) 

(9) Telescopio Carlos Sanchez (Canary I., Spain) 
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FIGURE CAPTIONS 


cloud (IC) from July 19 to September 23, 2009, in the visible continuum wavelengths ~ 
400-650 nm where it appears as a dark feature against the background Jupiter;iclMds. 




Dates indicated. 
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Figure 1. Series of cylindrical map projections showing the evolution of the impact 


■ 


Figure 2. Series of cylindrical map projections showing the evolution of the impact 
cloud (IC) from July 20 to October 29, 2009, in the near infrared spectral range 
covering the wavelengths of the methane absorption bands at 890 nm and 2.12-2.3 pm 
where it appears as a bright feature. Dates are indicated in each panel. 


Figure 3. Zonal expansion in System 111 longitude of the east and west edges (extremes) 
of the impact cloud (IC). main spot and crescent. Circles: measurements in the near 


infrared (2.12-2.3 pm) from July 19 to October 29, 2009; Crosses: measurements in the 


visual range u; 





PW and HST images from 19 July to 20 September 2009. 


Figure 4. Longitudinal drift of individual features within the IC tracked in the visual 

■tm 

S.yX-y ’*• 

range from IOPW images. A linear least-square fit to the motion of each feature is 




s Shown. Solid black circles correspond to a reference white spot (identified as W Spot 
placed at latitude -51.05°) not pertaining to IC that showed an oscillation in its longitude 


motion. 
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Figure 5. Longitudinal drift of individual features within the IC tracked in the near 
infrared (2.12-2.3 pm). A linear least-square fit to their motion is shown for each 
feature. 

Figure 6. Latitudinal drift from latitude -56.0° to -53.5° of a single cloud element within 
the IC, from July 23 to August 5. The inset shows the IC and the arrow marks the 
direction of the meridional motion of the cloud element. 
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Figure 7. Left: Zonal wind speeds of each individual element (dots) tracked within the 
IC using the images obtained in the visual range from I OP W (red) and HST (black), and 
in the near infrared (blue). Comparison is made with the background zonal wind profile 
measured in 2007 using the HST images (solid black line with error bars) and with 
individual features (circles) tracked in the same HST images as the IC. The impact point 
is represented by a large green circle. Right: All the IC individual cloud motions (dots) 
and their fit to a simple parabolic wind profile function (red curve) compared with the 

HST 2007 profile (continuous line). 
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Figure 








: east and west expansion of the IC and clump formation in its interior 


Ju 


££ -OK 


to Aug. 3 as observed in the methane absorption bands that sense the IC 


top .aerosol level. The formation and tracking of different spots within IC are marked by 
the lines and arrows. From top to bottom: July 23 (HST, 890 nm), July 25 (Keck, 2.12 
pm), July 27 (TNG, 2.3 pm), August 3 (TNG, 2.3 pm), August 3 (HST, 890 nm). 


Figure 9. Long-term evolution of the IC surface reflectivity in the near infrared 
methane absorption band at 2.2-2.3 pm (crosses) showing the exponential decay fitted 




function (solid line). The rectangle marks the reference calibration (see text) and the 
horizontal line marks the cloud background reflectivity. 


Figure 10. Polar map projections of the south polar region of Jupiter from HST images 
in the ultraviolet (A, B, C) and methane absorption band (D, E, F) showing the impact 
(dark in UV, bright in the methane band) and the planetary waves at north and sbuth. 
Dates and wavelengths: (A) July 23, 275 nm; (B) August 3, 225 nm; (Q Au^ust 8, 275 
nm; (D) July 23, 890 nm; (E) August 3, 890 nm; (F) August 8, 890 nm. I i " 


Figure 11. Snapshots of the temporal evolution of the 1C after 40 days resulting from a 
simulation of the IC spread by a purely advection model using 1 1 vertically detached 
layers. Each panel shows the vertically integrated amount of debris material with a 
color code that maximizes the structure. The last panel shows material which is 1 0 times 
less bright than the first one. Stripes in this figure appear from the limited number of 
vertical layers used to represent tM vertical wind shear. 


Figure 12. Comparison between EPIC numerical simulations and HST observations of 

' 

the IC. Left: maps of the potential vorticity (PV) at the altitude level corresponding to 

■lllk .jW 

the 180 rabar isobaric surface in the Jovian atmosphere resulting after 15 and 25 days 
following an initial impulse thermal anomaly disturbance of 500 Wm : . Right: HST 
images of the aerosol impact cloud after 15.1 and 20.5. days from the impact (see also 
Figure 1). 


Figure 13. (a) Total particle optical thickness from 10 mbar to 200 mbar as a function 
of time. Solid line is used for scenario A and dashed lines for scenarios B (upper line) 




and C (bottom line), (b) Particle density versus time for upper layer (n upper ) and bottom 
layer ( niower )• As in (a), solid line shows scenario A. Note that scenario B yields the 
bottom dashed line for n up per and the upper one for ni ower , just the opposite to scenario C. 


See text for more details. 
























I. #**«£• 































5000 



5040 5080 $120 5160 

Julian Date Reduced (+2450000.00) 

19 July 2009 (12:00 UT) * 2455032.00 



apnjiBuoi III UiajsAs 


Julian Date deduce 











«'* % 








July 


23 


25 


July 


July 


27 


August 3 

















Ofsa 




















t 

lip 

f pj 


111 * 





48 


We present a study of the long-term evolution of the cloud of aerosols produced in the 
atmosphere of Jupiter by the impact of an object on 19 July 2009. 

The impact cloud expanded zonally from ~ 5000 km to 225,000 km (29 October), 
remaining meridionally localized within a latitude band from 53.5°S to 61.5°S 
planetographic latitude. 

We find that the westward jet at 56.5°S latitude increases its eastward velocity with 
altitude with vertical wind shear of 1 ms' 1 per scale height above the tropopause. 


A pure advection of the aerosols by the winds and their shears and a nonlinear 
simulation (EPIC code) of the evolution of the potential vorticity field generated by a 
heat pulse, reproduce the observed global structure of the cloud, j|sfopf§|al depth 
decrease, and the dominant zonal dispersion of the aerosols. 









